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SUMMARY 


Aerodynamic  influence  coefficients  for  symmetrically  loaded  wings 
of  arbitrary  plan  form  in  subsonic  flow  are  derived  from  a  simple 
empirical  method  of  estimating  spanvise  lift  distributions.  The 
application  of  the  coefficients  to  an  aeroelastic  analysis  is 
discussed. 


INTRODUCTION 


In  the  aeroelastic  analysis  of  reference  1  the  spanvise  lift  dis¬ 
tribution  was  assumed  to  be  given  by  strip  theory  with  over-all  reduc¬ 
tion  and  rounding  at  the  tip  to  take  aerodynamic  induction  approxi¬ 
mately  into  account.  The  spanwise  lift  distribution  may  be  incorporated 
into  the  aeroelastic  analysis  more  accurately  by  using  suitable  aero¬ 
dynamic  influence  coefficients.  In  this  paper  a  method  is  presented 
for  determining  approximate  aerodynamic  influence  coefficients  for 
symmetric  loadings  from  the  method  of  estimating  lift  distributions 
presented  in  reference  2.  An  approximate  means  of  taking  the  effect 
of  the  fuselage  into  account  is  also  presented  in  the  present  paper; 
it  is  based  on  the  assumption  that  the  fuselage  carries  an  amount  of 
lift  which  depends  solely  on  the  angle  of  attack  of  the  wing  at  its 
root . 


Although  this  paper  is  intended  primarily  as  an  adjunct  to  refer¬ 
ence  1,  the  method  and  coefficients  presented  herein  may  be  used  for 
calculating  lift  distributions  apart  from  reference  1.  In  order  to 
facilitate  the  application  of  the  present  paper  to  such  calculations, 
all  required  definitions  and  numerical  factors  are  presented  herein. 

As  in  reference  1  matrix  notation  is  used  in  the  present  paper.  A 
short  summary  of  matrix  notation  and  algebra  and  a  discussion  and 
derivation  of  the  pesrtinent  integrating  matrices  are  given  in 
reference  1. 
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SYMBOLS 


b 

bf 


c 

c 

ce 

cs 

CZ 

CLa 

w 

f 

g 

M’M’fle] 

Nj. 

H 

s 

y 

y* 

a 

a 

n 

T)* 


span 

span  less  fuselage  width 

chord,  parallel  to  plane  of  symmetry 

average  chord 

average  of  root  chord  and  chord  in  plane  of  symmetry 

chord  in  plane  of  symmetry 

local  lift  coefficient 

wing  lift -curve  slope 

fuselage  width 

"ideal"  function  (see  fig.  1) 
factor  defined  by  equation  (ll) 


integrating  matrices 
first  row  of  K^_  matrix 

aerodynamic  -  influence  -coefficient  matrix 
total  wing  area 

lateral  ordinate  measured  from  plane  of  symmetry 

dimensionless  lateral  ordinate 

angle  of  attack 
average  angle  of  attack,  defined  by  equation  (2) 
lateral  ordinate  measured  from  wing  root 

dimensionless  lateral  ordinate 


(b72 
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* 


X  taper  ratio 

A  '  angle  of  sweepback  at  quarter -chord  line 


Subscripts: 

a 

b 

r 

Matrix  notation: 


additional 

basic 

root 


square  matrix 
diagonal  matrix 

column  matrix 
row  matrix 

unit  matrix 


CALCULATION  OF  INFLUENCE  COEFFICIENTS 
Wing  Alone 


In  reference  2  the  spanwise  lift  distribution  at  subsonic  speeds 
for  any  angle -of -attack  distribution  is  obtained  by  considering  sepa¬ 
rately  the  additional  lift  distribution  and  the  basic  lift  distribution. 
The  additional  lift  distribution  is  calculated  from  the  average  of  the 
chord  distribution  and  an  "ideal"  function  f  (given  in  reference  2  and 
reproduced  in  fig.  1  of  the  present  paper ),  which  has  been  obtained 
empirically  for  various  angles  of  sweep: 


(1) 
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The  average  angle  of  attack  a  is  obtained  from 

o  pb/2  nl 

a  =  —  /  ac  d y  =  /  a  -2.  dy*  (2) 

SJ0  ^0  c 

where  y*  is  the  dimensionless  lateral  ordinate  defined  by 

„  y 

V*  =  b72 

Similarly,  the  basic  lift  distribution  is  obtained  from  one-half 
of  the  product. of  the  chord  distribution  and  the  angle -of -attack 
distribution  for  zero  lift  (a  -a): 

fr)b  -  |  f(a  - »  (3) 

Any  sharp  corners  in  the  resulting  distribution  must  be  rounded  off. 

Equations  (l)  and  (3)  may  be  combined  and  written  in  matrix 
notation  as 

(4) 

where  |lj>  represents  a  column  of  lfs.  The  integral  required  in  cal¬ 
culating  a  may  be  evaluated  by  means  of  an  integrating  matrix.  The 
first  row  of  the  Kp  matrix  of  reference  1  LKlJx  serves  to  perform 
this  operation  and  round  off  the  basic  distribution  at  the  same  time. 
The  six-point  matrix  (shown  in  table  1  of  the  present 

paper)  is  based  on  spanwise  stations  located  at  0,  0.2,  0.4,  0.6 ,  0.8 , 
and  0.9  of  the  distance  over  which  the  integration  is  performed.  Con¬ 
sequently,  the  analysis  of  the  wing  alone  is  based  on  stations 
corresponding  to  y*  =  0,  0.2,  0.4,  0.6,  0.8,  and  0.9- 

In  order  to  fulfill  the  requirement  that  the  integral  of  the  basic 
lift  distribution  is  zero,  the  value  of  a  may  be  determined  from  the 
relation 
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or 


o 


where,  for  the  special  case  of  the  wing  alone. 


(6) 


The  value  of  g  is  presented  for  various  values  of  the  taper  ratio  in 

figure  2.  The  case  of  the  wing  alone  corresponds  to  —  =  0. 

b 


Equation  (5)  may  be  substituted  in  equation  (4)  and  both  sides 
multiplied  by  —  to  yield 


where  the  aerodynamic -influence -coefficient  matrix  Q  is  defined  by 


(8) 


The  K]_  matrix  is  a  square  matrix  consisting  of  rows  all  equal  to 


The  elements  of  ^ 


K- 


and 


LkiJi- 


are  given  in  table  1,  the  factor 


of  2  being  included  to  avoid  the  necessity  of  a  final  multiplication  by 

that  factor.  In  equations  (7)  and  (8)  the  chord  at  the  plane  of  sym¬ 
metry  could  be  used  instead  of  the  root  chord  cr  since  the  two  are 
identical  for  a  wing  without  fuselage. 


In  calculating  the  Q  matrix  the  f  values  are  first  read  for 
the  given  sweep  angle  from  figure  1  at  y*  values  of  0,  0.2,  0.4,  0.6, 
0.8,  and  0.9;  these  f  values  are  then  entered  in  a  diagonal  matrix  as 
shown  in  table  2  for  an  example  wing  with  30°  sweepback.  A  value  of  g 
is  then  read  from  figure  2  for  the  given  taper  ratio  (X  =  0.5  for  the 
example  wing)  and  multiplied  into  the  f  matrix.  The  resulting  matrix 


is  postmultiplied  by  the 


matrix  of  table  1. 


(See  table  2  for 
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t 

the  result  of  this  operation  in  the  case  of  the  example  wing. )  The 
matrix  is  then  added  to  the  product.  The  resulting  matrix  is 

c 

postmult Iplied  by  a  diagonal  matrix  of  the  —  values  at  the  previ¬ 
ew 

ously  cited  values  of  y*  in  order  to  obtain  the  desired  Q  matrix. 

(See  table  2.) 


Wing  and  Fuselage 

For  the  purpose  of  calculating  aeroelastic  effects  the  presence  of 
the  fuselage  may  be  taken  into  account  in  several  ways.  In  one  method 
the  wing  root  is  considered  to  be  at  the  plane  of  symmetry  and  the  part  of 
the  wing  inside  the  fuselage,  to  be  infinitely  stiff.  In  this  case  the 
coefficients  obtained  in  the  foregoing  section  may  be  used. 

Another  method  consists  of  calculating  the  aerodynamic  influence 
coefficients  for  stations  corresponding  to  tj*  (rather  than  y*)  values 
of  0,  0.2,  0.4,  0.6,  0.8,  and  0.9*  In  that  case  an  assumption  has  to  be 
made  concerning  the  amount  of  lift  carried  by  the  fuselage  and  the 
effect  of  this  lift  on  that  carried  by  the  wing. 

For  untwisted  wings  it  is  commonly  assumed  that  the  fuselage  car¬ 
ries  the  same  lift  as  the  part  of  the  wing  it  intercepts.  In  this 
analysis  the  same  assumption  is  made  for  a  twisted  wing,  and  the  part  of 
the  wing  intercepted  by  the  fuselage  is  assumed  to  be  at  the  angle  of 
attack  of  the  wing  root.  The  value  of  a  for  this  assumption  may  then 
be  obtained  from  the  requirement  that  the  basic  lift  distribution 
integrates  to  zero: 


MCA  TN  2092 


7 


and 


g  = 


(11) 


Values  of  g  are  presented  in  figure  2  as  a  function  of  the  fuselage- 
width  and  taper  ratios. 


With  the  value  of  a  given  by  equation  (9),  equation  (4)  takes 
the  form  of  equation  (7)  where  Q  is  now  defined  by 


H-  *  i&t!] 


(12) 


Comparison  of  equations  (9)  and  (12)  with  equations  (5)  and  (8),  respec- 


tively^  indicates  that  the 


matrix  is  obtained  by  adding  the 


quantity 

matrix. 


w 


cr 


to  every  element  of  the  first  column  of  the 


When  the  Q  matrix  of  equation  (12)  is  used  in  equation  (7),  the 
CCZ 

values  of  —  will  be  obtained  at  stations  corresponding  to  t)*  =  0 , 

0.2,  0.4,  0.6,  0.8,  and  0.9  for  the  given  values  of  the  angle  of ' attack 
at  those  stations.  The  relation 


y*  = 


may  facilitate  determination  of  the  stations  of  interest.  In  calcu¬ 
lating  the  Q  matrix  for  the  wing -fuselage  combination,  values  of  the 

f  function  and  the  —  ratio  are  obtained  at  those  stations.  The 

cr 

average  chord  ce,  which  is  calculated  from  equation  (9),  is  used  to  obtain 


the  term 


the  — 

2 


2  b7  cP  which  is  added  to  each  element  of  the  first  column  of 

matrix  of  table  1  to  yield  the  — 

2 


K- 


matrix. 


This  matrix  is  premultiplied  by  a  diagonal  matrix  obtained  by 
multiplying  the  f  matrix  by  the  value  of  g  taken  from  figure  2  for 
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the  given  configuration.  The  ff1]  matrix  is  added  to  the  resulting 
matrix  and  the  result  is  postmultiplied  by.  the  matrix  to  yield 

the  desired  Q  matrix. 

APPLICATION  OF  INFLUENCE  COEFFICIENTS 


The  aerodynamic  influence  coefficients  calculated  as  outlined  in 
the  foregoing  section  may  he  used  to  calculate  the  lift  distribution  of 
wings  with  arbitrary  sweep  and  taper  for  any  symmetrical  angle -of -attack 
distribution.  (For  discontinuous  distributions,  such  as  those  caused  by 
flap  deflections,  the  equivalent  5  values  of  reference  3  may  be  used 
in  order  to  obtain  the  rounding  of  the  lift  distribution  at  the  point  of 
discontinuity  called  for  in  reference  2.)  Once  the  coefficients  have 

been  calculated  for  the  plan  form  of  interest  (with  or  without  fuselage), 

ccz 

the  values  of  the  loading  coefficient  - —  at  stations  corresponding 

cr 

to  0,  0.2,  0  A,  0.6/  0.6,  and  0.9  of  the  distance  from  the  wing  root  to 
the  wing  tip  are  obtained  by  premultiplying  a  column  matrix  of  the 
values  of  the  angle  of  attack,  at  those  stations  by  the  influence - 
coefficient  matrix. 

When  the  influence  coefficients  are  used  in  conjunction  with  the 
method  of  aeroelastic  analysis  presented  in  reference  1,  steps  (6)  and 

(ll)  in  table  VI  of  reference  1  must  be  modified.  The  influence - 

coefficient  matrix  replaces  the  —  matrix  used  in  step  (l^  .  A  matrix 

cr 

must  be  calculated  (from  the  given  values  of  the  moment 
sub  c 

arm  ep  pertinent  to  subsonic  speeds  and  the  given  chord  ratio  — ^ 

and  postmultiplied  by  the  influence -coefficient  matrix.  The  resulting 
matrix  replaces  the  — — \  matrix  called  for  in  step  @  .  This 

LflrVcr/Jsub 

procedure  is  not  applicable  to  the  supersonic  case  since  the  aerodynamic 
influence  coefficients  have  been  obtained  only  for  the  subsonic  case. 
Elsewhere  in  the  method  of  reference  1,  the  lift -curve  slope  0L<x  of 

the  present  paper  is  used  in  place  of  the  effective  slopes  ]%  and 

me  I  ... 

nu  the  ratio  -  is  taken  as  1,  and  the  equivalent  angle  of  attack  a 

^®1*  __ 
of  reference  1  (not  to  be  confused  with  the  average  angle  of  attack  a 
of  the  present  paper)  is  merely  the  sum  of  the  geometric  and  structural 
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angles  of  attack,  A  convenient  way  of  obtaining  the  lift  coefficient 
for  any  calculated  loading  consists  of  multiplying  the  a  values 
obtained  from  equation  (5)  or  (9)  by  the  lift-curve  slope  Ci^.  In  all 
other  respects  the  calculations  indicated  in  reference  1  are  unaffected 
by  the  use  of  the  influence  coefficients. 

Compressibility  corrections  may  be  applied,  as  suggested  in  refer¬ 
ence  2,  by  using  an  f  function  corresponding  to  an  effective  sweep 
angle  and  by  using  a  corrected  lift-curve  slope. 


DISCUSSION 


The  influence -coefficient  matrices  of  this  paper  constitute  an 
expression  of  the  empirical  method  of  reference  2  in  matrix  form.  The 
matrix  elements  are  influence  coefficients  in  the  generalized  sense  that 
if  multiplied  by  the  angle-of-attack  values  along  the  span  and  summed  up 
they  will  yield  the  value  of  the  lift  at  a  given  station.  Individually 
the  elements  do  not  necessarily  represent  the  value  of  the  lift  at  a  given 
station  due  to  a  unit  angle  of  attack  at  another  station  in  the  more 
narrow  meaning  of  the  term  "influence  coefficient." 

The  discussion  in  reference  2  of  the  merits  and  demerits  of  the 
empirical  method  applies  to  the  aerodynamic  influence  coefficients  as 
well.  In  short,  at  subcritical  speeds  they  may  be  expected  to  yield 
results  which  are  more  reliable  than  those  of  a  strip-theory  analysis 
but  not  so  reliable  on  the  average  as  the  results  of  one  of  the  more 
refined  theories  for  spanwise  load  distribution.  The  results  differ 
from  those  of  the  method  of  reference  2  in  only  one  respect.  The 
rounding  of  the  basic  loading  curve  required  in  reference  2  is  performed 
by  integrating  matrices  in  the  present  paper.  In  order  to  show  the 
effect  of  this  difference,  the  results  of  the  two  approaches  for  the 
basic  lift  distribution  of  an  untapered  unswept  wing  with  washout  of 
1  radian  are  shown  in  figure  k.  Also  shown  are  results  calculated  by 
means  of  Multhopp's  and  WeissingerTs  analytical  methods  and  on  the  basis 
of  the  modified  strip  theory  used  in  reference  1. 

The  influence -coefficient  method  and  the  method  of  reference  2 
yield  almost  identical  results,  which  are  in  fair  agreement  with  the 
results  of  Multhopp’s  and  Weissinger’s  analytical  methods.  The  agreement 
of  the  results  calculated  by  the  modified  strip  theory  of  reference  1  with 
the  analytical  results  is  somewhat  poorer  on  the  average.  (The  results  of 
uijmodified  strip  theory  would  be  much  poorer;  for  the  case  of  fig.  h  the 
loading- coefficient  curve  would  be  80  percent  larger  at  every  point  than 
that  of  the  modified  strip  theory  and  would  not  be  rounded  off  at  the  tip.) 


The  methods  used  for  the  calculations  represented  in  figure  k 
(except  that  of  modified  strip  theory)  would  yield  almost  equal 
results  for  the  additional  loading,  which  would  he  intermediate  between 
elliptical  and  rectangular;  whereas  the  result  of  modified  strip  theory 
would  be  a  rectangular  distribution  rounded  only  at  the  tip.  The 
results  of  the  influence -coefficient  method  are  therefore  superior  to 
those  of  the  modified  strip  theory  for  both  basic  and  additional  lift 
distributions,  although  the  difference  is  not  so  large  for  tapered  wings 
as  for  the  untapered' wing  considered  in  figure  4.  Since  the  agreement 
of  the  results  of  the  influence -coefficient  method  with  those  of  the 
analytical  methods  is  somewhat  poorer  for  basic  lift  distributions  than 
for  additional  lift  distributions,  the  influence  coefficients  may 
be  expected  to  be  most  useful  when  the  angle  of  twist  is  less  than  the 
angle  of  attack.  The  influence  coefficients  will  probably  furnish 
fairly  reliable  results  for  the  change  in  loading  due  to  aeroelastic 
effects  but ' somewhat  less  reliable  results  for  the  divergence  speed. 

In  general,  they  tend  to  overestimate  the  lift,  particularly  at  the  wing 
tip,  so  that  aeroelastic  effects  are  overestimated  to  some  extent.  For 
similar  reasons  the  influence  coefficients  of  this  paper  are  not  appli¬ 
cable  to  antisymmetric  lift  distributions  because  these  lift  distri¬ 
butions  do  not  have  an  additional  part. 


CONCLUDING  REMARKS 


A  method  has  been  presented  for  calculating  aerodynamic  influence 
coefficients  for  symmetrical  loadings  at  subsonic  speeds  from  a  simple 
empirical  method  of  estimating  spanwise  lift  distributions.  The  coeffi¬ 
cients  are  particularly  suited  for  aeroelastic  analyses  and  may  be 
expected  to  give  good  results  for  the  changes  in  spanwise  loading 
associated  with  aeroelastic  effects  but  probably  less  reliable  results 
for  the  divergence  speed. 
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TABLE  1 

VALUES  OF  THE  I  K-,  I  .  ifel  ,  AND  %Txl  MATRICES 


0.5000 

0 

0 

0 

0 

0 

0 

.5000 

0 

0 

0 

0 

0 

0 

.5000 

0 

0 

0 

0 

0 

0 

.5000 

0 

0 

0 

0 

0 

0 

.5000 

0 

0 

0 

0 

0 

0 

.5000 

Factor 


Loading  coefficie 
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5^ 


Dimensionless  lateral  ordinate,  y* 

.  Figure  4.-  Basic  lift  distribution  of  rectangular  wing  of  aspect 
ratio  5  with  unit  linear  symmetric  twist  in  incompressible 
flow. 
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